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We present a theoretical study of transport properties through superconducting contacts based 
on a new formulation of boundary conditions that mimics interfaces for the quasiclassical theory of 
superconductivity. These boundary conditions are based on a description of an interface in terms of 
a simple Hamiltonian. We show how this Hamiltonian description is incorporated into quasiclassical 
theory via a T-matrix equation by integrating out irrelevant energy scales right at the onset. The 
resulting boundary conditions reproduce results obtained by conventional quasiclassical boundary 
conditions, or by boundary conditions based on the scattering approach. This formalism is well 
suited for the analysis of magnetically active interfaces as well as for calculating time-dependent 
properties such as the current-voltage characteristics or as current fluctuations in junctions with 
arbitrary transmission and bias voltage. This approach is illustrated with the calculation of Joseph- 
son currents through a variety of superconducting junctions ranging from conventional to d-wave 
superconductors, and to the analysis of supercurrent through a ferromagnetic nanoparticle. The cal- 
culation of the current- voltage characteristics and of noise is applied to the case of a contact between 
two d-wave superconductors. In particular, we discuss the use of shot noise for the measurement of 
charge transferred in a multiple Andreev reflection in d-wave superconductors. 



I. INTRODUCTION 

At the end of the 1960's several authors demonstrated 
that the complete standard theory of superconductivity, 
both in equilibrium and out of equilibrium, can be for- 
mulated in terms of a quasiclassical transport equation!!]. 
This theory, known as quasiclassical theory of super- 
conductivity, combines Landau's semiclassical transport 
equations for quasiparticles with the concept of pairing 
and particle-hole coherence that are the basis of the BCS 
theory. The quasiclassical theory provides a full descrip- 
tion of superconducting phenomena ranging from inho- 
mogeneous superconductors to superconducting phenom- 
ena far from equilibrium. In its traditional form, the 
quasiclassical theory of superconductivity is restricted 
to the description of weak perturbations of a supercon- 
ductor: the external perturbations (magnetic field, vari- 
ations in the chemical potential, etc) should be weak 
(V <C Ep), where Ep is the Fermi energy, of long wave- 
length (q kp 1 ), where kp 1 is the Fermi wavelength, 
and of low frequency {fko <C Epp. Although this limit 
covers many situations of interest in superconductivity, 
some important cases fall outside the range of validity of 
the usual scheme of the quasiclassical theory: strong im- 
purities, walls, interfaces, etc. In the case of the analysis 
of the transport properties of superconducting contacts, 
the quasiclassical theory is complemented by appropri- 
ate boundary conditions, the so-called Zaitsev boundary 
conditional, or by its generalization to magnetically ac- 
tive interfaces by Millis, Raincr and Saulsu. Although 



these boundary conditions provide a formal solution of 
the problem of a strong perturbation due to interfaces, 
their highly non-linear form introduces many problems, 
e.g. these boundary conditions have spurious solutions 
which require special care, in particular in numerical im- 
plementations. At this point, we should mention that the 
recent progress in the solution of these boundary condi- 
tions make them much more tractableoQ. 

In the field of electronic transport through supercon- 
ducting junctions there are many basic situations in 
which a complete understanding is still lacking. This is 
partially due to the difficulty of applying the quasiclassi- 
cal theory because of the lack of simple and manageable 
boundary conditions. Most of these situations are re- 
lated to either the current-voltage (I-V) characteristics 
of contacts between two superconductors with arbitrary 
transmission, or to situations in which there are magneti- 
cally active interfaces giving rise to spin-dependent trans- 
port. For instance, in the context of high temperature 
superconductors (HTS),r-svhich by now are believed to be 
d-wave superconductors^!, we can mention N-I-S tunnel 
junctionsLTtj, grain boundary junctions, widely used for 
the realization of Joaephson junctional3il3, and YBCO 
S-N-S edge junctionsEJjlLM In the two latter systems the 
description of the I-V characteristics reouires an analysis 
of multiple Andreev reflections (MAR)Ej, while all sys- 
tems above require a careful self-consistent study of the 
superconducting state in the presence of pair-breaking, 
caused by bulk or surface impurities and by the sur- 
face itself, in order to account for results obtained in 
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experimentJHll by theory&El. Other examples, which 
take place in the context of conventional superconduc- 
tors, are the S-N-S point contacts, where the normal re- 
gion N has a length comparable to the superconducting 
coherent lengthMEa. These systems demand from the 
theory a detailed analysis of the interplay between the 
proximity effect superconductivity in the normal metal 
and the occurrence of MAR. This analysis is currently 
in progress making use of the approach described in this 
worlH. 

With respect to the transport through spin active in- 
terfaces there have been an extensive experimental ef- 
fort to -expJorc tunneling through magnetic insulator 
barriersr 7 ^ 28 ! by probing the proximity effect in S/E=. 
structures£j'E!l, and by constructing S/F-multilayersEil 
(see also references therein). Andreev reflection us- 
ing ferromagnet/superconductor point contacts have 
also been used-, to probe the spin polarization of the 
ferromagnetoH Finally, supercurrents have been re- 
ported in S/F/S_iunctions by Veretennikov et no and 
by Gandit et a£j. To model an S/F/S-junction one 
can follow one of two routes. The first route is to as- 
sume an extension of a ferromagnetic metal, now char- 
acterized by a length and an exchange field, separat- 
ing the two superconductors. pEollowing this route, 
both critical current oscillationsE§EZl and the effect, p£ 
the exchange field on the Andreev bound statesc§E2l 
have been studied. The limitation here is that the ap- 
proach is restricted to small exchange fields in the fer- 
romagnet, i.e. fields that are comparable to the su- 
perconducting gap. The second route, which we sub- 
scribe to, is to treat the effect of the ferromagnet as a 
boundary pxabLep. connecting the two superconducting 
half-spacesEH'QEJ. Using Bogoliubov-deGennes equa- 
tions and a WKB-approach for the ferromagnetic bar- 
rier, i.e. the scattering approaches, Josephson .cuxrent- 
phase relationscij and quasiparticle tunnelingac3 have 
been studied for both conventional s-wave and uncon- 
ventional d-wave superconductors. 

Parallel to the quasiclassical treatment of interfaces, 
a different formalism referred to as the Hamiltonian ap- 
proach has emerged as a new interesting tool for the anal- 
ysis of transport through superconducting contactsOuJ. 
This approach is based on modeling the contact by a 
simple Hamiltonian in combination with non-equilibrium 
Green's function techniques^. The origin of this ap- 
proach can be traced back to the work of Bardeen who 
introduced the tunnel Hamiltoniafl-,approximation in or- 
der to describe a tunnel junctionea. Most of the calcu- 
lations of current through superconducting contacts in 
the early 60's were based on this tunnel Hamiltonian, 
and were restricted to the lowest order transport pro- 
cesses such as the calculation of the Josephson current 
in a S-I-S junctionE3i£j. Multi-particle— tunneling was 
first discussed by Schrieffer and WilkinsLHl in their mul- 
tiparticle tunneling theory as a possible explanation for 
the observed sub-gap structure in superconducting tun- 
nel junctions. The contributions of these higher order 



processes were found to be divergent, which has led to 
the quite extended belief that the Hamiltonian approach 
is pathological except for describing the lowest order tun- 
neling processes. Howevec, in a series of works, Caroli 
and coworkers in the 70's^.and, more recently, Martfn- 
Rodero and collaboratorsCJO have shown that one can 
get rid of all the old pathologies of this approach by using 
a local representation and by summing the series of tun- 
neling processes up to infinite order. In particular, within 
this technique a microscopic theory of Multiple Andreev 
reflections has been developedc3. This theory describes 
quantitatively the I-Y characteristic of superconducting 
atomic-size contaciO'Ei 

In this paper we show how the Hamiltonian ap- 
proach may be brought into the quasiclassical theory 
by integrating out large energies right at the onset. 
The resulting boundary condition reproduces results ob- 
tained by thej-conventional quasiclassical Zaitsev bound- 
ary conditionEra or by, boundary conditions based on the 
scattering approachE3 t2l. The boundary condition is 
also well suited for calculating time-dependent properties 
such as the current-voltage characteristics or the current 
fluctuations of S-I-S junctions. The paper is organized 
in the following manner: In section II we give the energy 
integration of the Hamiltonian approach and state the re- 
sulting boundary condition. In section III we show how 
the current through a contact may be calculated from the 
boundary t-matrix. In this section we also calculate the 
Josephson current resolved in energy and on trajectory 
for different types of superconductors and for different 
types of coupling between the two superconductors. In 
section IV, we discuss the boundary condition at a fi- 
nite bias applied between the superconductors. This is 
then applied to the case of two coupled d-wave supercon- 
ductors. Finally, in section V, current-fluctuations are 
discussed and the current-current correlator is derived 
within our method. This is then applied to compute the 
trajectory resolved current-fluctuations of two coupled d- 
wave superconductors. 



II. DESCRIPTION OF THE APPROACH 



The system of study is two semi-infinite superconduct- 
ing electrodes coupled over some type of interface bar- 
rier. Our approach is to artificially separate the problem 
in to two parts in order to pose a boundary condition for 
the interface. The first part consists of calculating the 
Green's function of either conductor, extending to ±00 
respectively, in presence of a hard surface at x =j-iO. To 
this part of the problem, the quasiclassical theorycl is our 
theory of choice. It has been shown that strong perturba- 
tions, such as rigid walls may be included into quasiclas- 
sical theory by means of effective boundary p^nnditions 
posed for the quasiclassical Green's functiorBErEj. To 
couple the two electrodes, from now denoted left (L) and 
right (R) , we assume a phenomenological Hamiltonian as 
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follows 
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The potentials ul.r and url, with ujj L = vlr = V, 
act as hopping elements connecting the two electrodes 
L and R. The perturbation given by Ht is short ranged 
(~ Af <C Co) and may be strong (v ~ The local 

character of Ht allows us to view it as a single strong 
impurity in case of a point contact or as a line of strong 
impurities in case of an extended contact between the 
two electtpdes. Following the work of Thuneberg and co- 
workera23o the single impurity or, in case of the line,-e|f 
impurities following the work of Buchholtz and RainerEil, 
this strong perturbation may also be incorporated into 
quasiclassical theory via a T-matrix equation. Anticipat- 
ing the result for the T-matrix, the effect of the contact 
between the two electrodes on the physical quasiclassical 
Keldysh-Nambu matrix Green's function, or propagator, 
<7i, in electrode i = (L, R) enters as a source term in the 
transport equation for £/j(j£,) along a trajectory q, 



which together with the Dyson equation 



(2) 



Here R c is the position of the contact and tfc is the Fermi 
velocity at point q, on the Fermi surface. The Green's 
function <7oo,i is an intermediate Green's function ob- 
tained by solving the hard wall boundary condition of 
the separate electrodes, i.e. without taking the contact 
into account but using the self energies q and A, eval- 
uated using the physical propagator, cji, satisfying Eq. 

Our objective is to find the quasiclassical T-matrix, 
t, giving the source term in the transport equation (j^) 
above. The starting point is a conventional many-body 
perturbation theory for the Hamiltonian Ht- To proceed 
we artificially enlarge our Hilbert space with a " reservoir 
quantum number" (L, R) and the functions entering are 
the matrices 
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The matrix elements are the usual Keldysh-Nambu ma- 
trices of non-equilibrium superconductivity!! Especially, 
the Green's functions Gqo.l and Goo,i? are the Green's 
functions for the uncoupled left and right electrode. The 
coupling elements, vlr,rl, between L and R are propor- 
tional to the unit matrix in the Keldysh space and in 
Nambu space adopt the form 
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constitutes a closed set of equations that are to be 
brought into quasiclassical form. We have given two dif- 
ferent ways of summing the series which correspond ei- 
ther to "dressing" the perturbation (equations (||) and 
(||)) or to "dressing" the Green's function (equations (^) 
and (0)). The two sets of equations are equivalent and 
two useful relations 
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With this, we write the T-matrix equation 



follow directly. Here, and above, the o-product is short- 
hand for integration or summation over common argu- 
ments. Starting from equation (Q), using equation (|^) 
and the second of the two relations (0) , it is straight for- 
ward to get the following closed set of equations, closed 
separately fop-pne and each of the components, Ty, of 
the T-matrixEl, 

Tll = v L R ° Goo, r ° v RL 

+ vlr ° Goo,i?, ° vrl o Goo,l ° T L l 
Trr = vrl o Goo.L ° v L r (8) 

+ VRL ° Goo.L O V L R O Goo,_R O Trr 

Tlr = v L R + v L r ° G^r o v rl ° Goo.l ° T L r 

Trl = Vrl + Vrl o Goo,L O V L R o Goo.i?, o f RL . 

The equations above depend only on the Green's func- 
tions Goo,l and Goo.fl of the two uncoupled systems. 
Since the full Green's function Gij has been eliminated 
from the T-matrix equations there are no Green's func- 
tions with spatial arguments in both systems. Together 
with the short range of vlr,rl this means that we 
can directly perform the quasiclassical ^-integration on 
the T-matrix equations and substitute the quasiclassical 
Green's functions, <7oo,i, for the full ones, Goo,i, above. 

At the quasiclassical level, the Green's functions 
goo,i(i^.]t,t') at the interface in equations (||) depend on 
the position on the Fermi surface, fy, and of two times 
(t, t'). The coupling elements will be assumed to be time 
independent but may couple different points and f(, on 
the Fermi surfaces of the two conductors. The exact form 
of the (/£,, ^-dependence of vlr is a degree of freedom 
in the model that allows us to consider different types of 
transport through the interface. We can now write down 
the equations for the quasiclassical i-matrix components 
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tLL = {vLR®9oo,R<SVRL)j£ 

+ ((vLR»goo,R®VRL®goo,L®tLL}f£')f£'' 
tRR = (#R£®5oo,£8>Ukfl)f£' (9) 

+ \\VRL®9co,L®fiLR®9oo,R®'bRR)ig)ty!' 
thR = VLR + ((vLR®goo,R®VnL®goo,L®tLR)p£')j^'' 
tRL = Vrl + ((vRL®goD,L®VLR®goo,R®tRL}f^'}j^", 

where we suppressed the explicit dependence of the func- 
tions on jj, and on time variables. The earlier o-product 
in Eq. (||) is replaced by the (^-product in the quasiclassi- 
cal expression. The O-product stands for an integration 
over a common time variable together with a normal ma- 
trix multiplication in the combined Keldysh-Nambu and 
spin spaceB. A left-over from the ^-integration is the in- 
termediate averaging over position on the Fermi surface 
as indicated by (• • >)fe — J FS ■ ■ ■ d%. The quasiclassi- 
cal i-matrix entering in to equation (0) is the forward 
scattering limito 

tij tij {ik libi^i^ ) 

of Eq. (|jj) and in general it depends on two times (t,t'). 

Following Ref. @, let us summarize the procedure for 
calculating the quasiclassical propagators in the presence 
of an interface: 

1 . To find g x , we solve the conventional quasiclassical 
equations, the Eilenberger equation or the Usadel 
equation, for the uncoupled electrodes in equilib- 
rium using hard- wall boundary conditions. 

2. Use (?oo to solve the quasiclassical T-matrix equa- 
tions (H). 

3. Solve the inhomogeneous quasiclassical equation 
(||) for the physical propagator g. 

4. Use g to calculate the "smooth" self-energies li and 
Aj which enter the quasiclassical equations for g^ 
and for g. 



The whole scheme amounts to a set of linear differen- 
tial equations for g^ and g, coupled in a non-linear way 
by the T-matrix and the self-energy equations. Substan- 
tial simplifications can be achieved in the case of low 
transmissive tunnel barriers or point contacts. In these 
cases one can neglect the influence of the neighboring 
electrodes in the calculation of the self-energies and then 
the equations for goo and ii and A^ decouple. 



III. JOSEPHSON CURRENTS 

As a first application of the boundary condition we 
calculate supercurrent through a variety of contacts con- 
necting two superconducting reservoirs. The current con- 
tribution from a given trajectory, fy, and at a given 
energy, e, may be calculated directly by integrating 
the transport equation (||) along the direction given by 
■^(j^,). This is easily seen as on the considered trajectory 
away from the contact the physical propagator giipp) co- 
incides with the intermediate propagator <7oo,i(f^) calcu- 
lated by the impenetrable surface condition. In absolute 
vicinity of the contact only the source term in Eq. (|^), 
[ia(fy,,i±,), goo : i(f^,)]S(R — R c ), contributes and results in 
a jump in the Green's function. The magnitude of this 
jump is given by integrating Eq. (|^) over the interval 
]0_, 0+[. Performing the integral results in the scattered 
propagator 
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9i+(n) = &-(&) — (!0) 

V F COS (pi 

where (pi is the angle tj?(ftO makes with the contact nor- 
mal. Note that jji-fy) = 3oo,i(ft) and to calculate 
i- c - t ne propagator along the trajectory (f^,) cou- 
pled to (fy) by pure surface scattering we must solve for 
gi+iii) along the path given by Urifj). The propagator 
at the contact, computed in (|10|), may now be inserted 
in the current formula 



3(T) = eNp J TrMky%; £ % = eN F j A (11) 

with Np the density of states at the Fermi level in the normal state. Since the Josephson current is an equilibrium 
property, the Keldysh Green's function components of g are in this case simply related to the Retarded (R) and 
Advanced (A) ones as g K = (g R — g A ) tanh(e/2T), and with (g A (jj.; e)y — f 3 g R (j^,; e)f 3 , we get the energy and 
trajectory resolved current contribution across the contact, evaluated in the left superconductor, as 

tanh(^). (12) 

The lonely first intermediate Green's function, 5oo,i(ft-), in Eq. ( fTo| ) explicitly drop out of the current in the angle 
average since they obey the impenetrable surface boundary condition. 

So far no reference have been made to the modeling of the contact and the t-matrix element t R L . The contact 
model depends on the choice of the momentum dependence of the coupling elements, vr RL = vlr. — v{j^.,j^.). Two 
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extreme models for the (f^,, ^-dependence will be considered: a totally disordered contact, v(f^,f^.) — v, i.e. the 
coupling across the contact retains no memory of the momentum direction, and a momentum conserving contact with 
ti(g,,g!) = v8{% — The i-matrix equations for the two types of contact, dropping superfluous indexing, read 



(13) 



for the disordered contact and 



(14) 



for the momentum conserving contact. For either model, the i-matrix equation above is simple to invert after inserting 
the retarded Green's functions gn(%) and 
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with the phase difference x across the junction and the upper (lower) signs of the phase refer to the right (left) 
electrode. In equilibrium the T-matrix equation is simply an algebraic equation in energy space, which can be 
trivially inverted. The energy and trajectory resolved current is written 
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tanh( 
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(16) 



Here we have traded in the coupling strength, v, for the 
transmission coefficient T>. The two are simply related 
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j(T) = eN F V- ; — T tanh( — — — ). (19) 



[1-X>sin 2 f] 



2T 



In equation (|l^) it should be noted that a second tem- 
perature dependence enters via the the temperature de- 
pendent gap, A(T). 



The scalar coupling element, v, is from now on a di- 
mensionless quantity. We have normalized the original 
coupling element in units of (1/ttNf) to get rid of the 
different prefactors that appear in the angle averages in 
the T-matrix equation ([)]). 



A. Josephson current between two s-wave 
superconductors 

Assuming that the two electrodes both are s-wave su- 
perconductors, we have the Green's functions on either 
side of the contact 



g{%)R(L) = 



e Ae ±l */ 2 
-Ae^*/ 2 -e 



(18) 



where = [A 2 — e 2 ] i . Since the s-wave superconductor is 
isotropic it does not matter which model for the coupling, 
( |l3| ) or (|l4|), we choose. Additionally, in the absence of 
surface depairing effects it is sufficient to know the bulk 
Green's functions ( |l8| ) to calculate the current contribu- 
tion (12). Using Eq. ( |l6| ) we find the known result JJaal 
the Josephson current is carried by junction statesQEHl 
located at £j(x, T) = ±A(T)[1 - X>sin 2 (x/2)]i The 
total current is the sum of all contributions and reads 



B. Josephson current between two d-wave 
superconductors 

To emphasize the modeling of the (j^,, ^-dependence 
of the coupling across the junction and the importance 
of using the correct surface Green's functions, we now 
study the current-phase relation of two d-wave supercon- 
ductors. A realization of a d-wave order parameter is 
A^_, = Acos2(0^_, — a). The magnitude and sign of A^ 
depends on the position on the Fermi circle and this is 
measured by the angle, 4>j^ j the angle % makes with the 
crystal a-axis. The angle a tracks the relative junction- 
to-crystal a-axis orientation. If a = ±ir/4 and specular 
quasiparticle scattering at the interface is assumed the 
order parameter seen along a trajectory changes sign at 
the surface and an Andreev bound state forms at zero en- 
ergy for every trajectory jj, [ We will stick with the 
junction realization a = ±7r/4. To incorporate the effect 
of these surface states into the current-phase relation one 
must use the surface Green's functionsEj'tll 



where fift = [A? — e 2 l s . The factor so discriminates 
between two types of junctioncSl'ESl: sr = —1 (an = at) 
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is referred to as a symmetric, and sr = 1 (or — as 
a mirror junction. The convention for signs of the phase, 
X, are as for the s-wave superconductor. It should be 
said before proceeding that we are neglecting the pair- 
breaking effect of the surfaced and we assume constant 
order parameters up to the interface. This is for the sake 
of simple illustration and for the comparison of analytical 
results with other boundary conditions. 

Turning to the i-matrix equations and starting with 
the diffusive model of the point contact, one immediately 
find that the Josephson current is zero. This follows from 
the vanishing average, (A^)^ = 0. Due to this prop- 
erty the anomalous propagators vanish and therefore the 
Josephson current as well (see Eq. (Q). In the oppo- 
site limit, using the momentum-conserving model (14), 
the Josephson current is not zero. After inverting the t- 
matrix equation and evaluating the commutator in (12), 
we find the energy resolved current at f^. 



je(ih) = ± Im 



8m X 



£j(x;»0 



tanh(— ). (21) 



The sign of j s {p^) is (+) f° r a mirror and (-) for a sym- 
metric junction. As in the case of the s-wave junction we 
have junction states carrying the Josephson current. The 
position of these states depend on the type of junction in 
the following way 



£j( X ;fr) = ±Vv\A : 
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symmetric. 



(22) 



These j-utiction states were first found by Riedel and 
BagwellEJ using a_ scattering approach and, indepen- 
dently, by BarashE^ -using what is known as the Zaitsev 
boundary conditions. Performing the integral over the 
energy we write the trajectory resolved current-phase re- 
lations 



( cos \ tanh( — S '" 2 ) 
j(B.)=±27n/D^|x{ ' _ (23) 



sin ¥ tanh(^— J-fe! — — ) 



for the two junction types, the upper being the mirror 
and the lower the symmetric one. The total current is 
the trajectory average of j(f^) multiplied by eNp. 

The main purpose of this section was to show that the 
quasiclassical version of the point contact coupling of two 
electrodes is simple to use and can recover results known 
in literature. In case of unconventional superconduct- 
ing electrodes results depend in a crucial way on how the 
contact is modeled. This gives the Hamiltonian boundary 
condition an advantage in flexibility to the conventional 
Zaitsev boundary condition which coincides with the mo- 
mentum conserving contact (J14|) introduced above. 



C. Josephson current through a spin active interface 

As another illustration of the flexibility of this method 
we extend the discussion to currents through spin active 
interfaces, i.e. interfaces which flip the spin of the inci- 
dent electrons either by spin-dependent scattering within 
the interface, or by a difference in spin-orbit coupling on 
either side of the interface. The general boundary con- 
ditions that connect the quasiclassical propagators for 
superconducting metals across magnetically activa-inter- 
faces were introduced by Millis, Rainer and SaulsQ. Re- 
cently one of the authors^ derived the explicit solution 
of these boundary conditions for equilibrium Green func- 
tions. In order to compare with this solution, as in Ref. 
[ we shall analyze in this section the Josephson cur- 
rent through a contact of two isotropic s-wave supercon- 
ductors connected through a small magnetically active 
junction. 

In order to accommodate the spin dependence we en- 
large our space in such a way that every quantity is now 
a 2 x 2 matrix in spin space. In particular, the coupling 
elements are spin dependent and adopt the following gen- 
eral form 
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v 
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where v — ( W ^ ^ 

W 1,T v l,l 



(24) 



Let us stick to the case of S/F/S junction analyzed in 
Ref. [ Qj. In this case F stands for a small ferromag- 
netic particle or grain. This ferromagnetic material is 
treated as a partially transparent barrier which trans- 
mits the two spin projections differently. For spin-active 
interfaces the different components of the S'-matrix, Sy 
are 2x2 spin matrices. To proceed further a specific 
S'-matrix was chosen in Ref. [ |tJ to model the magnetic 
barrier 



S 



f S n 


S12 \ 






{ S 21 


S22 J 


= 





exp(ieer 3 ), 



(25) 



where Oj notes the Pauli-matrices spanning spin space 
and parameters (t, r) are the usual transmission and re- 
flection coefficients. The S-matrix (p5|) is one of the sim- 
plest choices that allows a variable degree of spin mixing 
at the interface and the spin mixing is parameterized by 
the spin- mixing angle 0. By this construction S only vio- 
lates spin conservation, i.e. it does not commute with the 
quasiparticle spin operator a. The angle will be con- 
sidered as a phenomenological parameter independent of 
the trajectory direction (for more details see Tokuyasu 
et al Ref. [ 71 1). Within the approach presented of this 
paper, one can easily model the previous S-matrix with 
a spin dependent coupling 
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exp(i6(73) 

exp(— 46(73) 



(26) 
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Again the quasiclassical surface Green functions are the inputs in this approach, i.e. the Green's functions at the 
interface calculated with impenetrable wall .boundary conditions. The spin-active boundary condition must be used 
at the contact also for the reflecting surfaceQ. In this simple case of a magnetic barrier the resulting spin-dependent 
propagators can be written in a 4 x 4 block-diagonal form 



( 5TT hi \ 

fn hi 

gi L -/; T 

V ~f n .g TT j 



(27) 



where the electron and anomalous parts of 0) in the upper left corner of guock{%] 0) can be written as 

e fl cos(9/2) -fisin(0/2) 



<m(©) 
/u(©) 



V* sin(0/2) + ficos(0/2) 

A e -'( e= Fx)/2 

£ fi sin(0/2) + ftcos(0/2) 



(28) 



/it( ) %H s in(0/2)+ficos(0/2) 

for trajectories with (f^, ■ n > 0). For trajectories with reversed momentum, i.e. (jj, • n < 0), the phase factor 
cxp[=Fi(6 =p x)/2] for functions fn(Q) and /xt(0) goes to exp[±i(0 ± x)/2]. Above, as in the earlier examples, 
the phase difference, x, between the two reservoirs is included and the upper (lower) signs refer two the right (left) 
reservoir. The components of the propagator in the lower right corner of guock are simply related to those in the 
upper left corner by the replacement Q — > ~ 0. 

The angle induces a mixing of the two otherwise separated spin bands. It is easy to see that the resulting density 
of states at the interface has Andreev bound states inside the gap. These states are located at £b,i(i) = ±Acos(0/2), 
with +(— ) for the spin-up (-down) branch. The existence of the sub-gap states alters the Josephson current-phase 
relation radically. The contribution to the current from the energy e, the trajectory y^, and spin band f (!) may 
be calculated directly from expression ([!(]). Keeping in mind that the phase that enters in Eq. ( |l6| ) is the phase 
difference over the contact wc write 



smx 



[CI cos(0/2) + esin(0/2)f - VA 2 sin 2 (x/2) 



(29) 



for the current carried by the spin-up band. The current carried by the spin-down band is given simply as 
ie,i(???! 0>x) = Je,T (ft-! — 0> X)- ^ a finite superconducting phase difference, the original two interface Andreev 
bound states are split up into four current carrying states located inside the gap at positions 



±A 



cos 2 (0/2) - X>cos(0) sin 2 (x/2) ± VX>sin(0) sin(x/2Wl - Psin 2 (x/2) 



1/2 



(30) 



These states give the total contribution to the current 
and their positions change from the tunnel regime, ej = 
±Acos(0/2), to ej = ±Acos[(x± ©)/2)] at perfect 
transmission. 
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FIG. 1. Zero-temperature supercurrent-phase relation for 
the S/F/S contact considered in this section. The four dif- 
ferent panels correspond to T> = 0.1,0.4,0.8, and 1.0. The 
spin-mixing angle Q is varied from top to bottom from to 
7r in steps of yjj . The dashed lines indicate the value of O for 
which the contacts become -k junctions. The supercurrent is 
normalized in units of the critical current density, jc, for the 
corresponding transmission and zero spin-mixing angle. 



FIG. 2. Energy resolved spectral current for the S/F/S con- 
tact with T> — 0.4 and 6 = 0.77T (see Fig. [l]right upper panel). 
This figure shows the total spectral current, sum of both spin 
contributions, as a function of energy according to Eq. (p9[), 
without the thermal factor. The different curves show the 
evolution of the four current-carrying Andreev bound states 
inside the gap with the superconducting phase difference. The 
curves are shifted vertically for clarity. 



In figure |J we show the current-phase relation for a 
set of transparencies V = 0.1,0.4,0.8, and 1.0. In each 
panel the spin-mixing angle is varied from to tt in 
steps of jq ■ As seen, for each value of T> there is a range 
> O c where the junctions are 7r-junctions. For small D, 
the critical spin-mixing angle O c is close to ir and with 
increasing T>, © c increases towards 5. The magnitude 
of the critical current is for all but the perfect trans- 
mission junction very asymmetric for and 7r junctions. 
In figure we show the energy resolved spectral current 
( p9| ) for V — 0.4 and 9 = 0.77T at different phase differ- 
ences, x, over the junction. At small phase differences 
the junction state initially at £j(x = 0) = ±Acos(0/2) 
splits in to two states carrying current in opposite direc- 
tions. This gives a small but positive current as seen in 
the corresponding current-phase relation in figure [l]. As 
the phase difference is increased the two states dispers- 
ing with phase towards e — from either side eventu- 
ally cross at e = 0. Above this phase difference, both 
current-carrying states at e < (e > 0) give current 
in the same direction and the magnitude of the current 
increases abruptly. 



To conclude this section, let us stress that these results 
reproduce the result obtained in Ref. ^, showing again the 
versatility of the boundary conditions introduced in this 
work. 
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IV. SOLVING THE T-MATRIX EQUATION AT 
AN APPLIED VOLTAGE 

In this section we shall describe how the T-matrix 
equation (^) can be solved in the case of a voltage-biased 
superconducting contact. As a constant bias V is applied 
across a junction between two superconductors, the phase 
difference oscillates with time according to the Joseph- 
son relation: <j>(t) — cfo + u>ot, where loq = 2eV/h is 
the Josephson frequency. This makes that every Green's 
function and T-matrix component depend on two time 
arguments. We show in this section how the time convo- 
lutions in the T-matrix equation can be handled. 

As we shall show in next section the current can be 
expressed, for instance, only in terms of the advanced 
and retarded components of thR- Thus, we concentrate 
ourselves in these components whose equations can be 
written as (see Eq. 9) 

^Zflt*) = ®LR + J dti J dt 2 VLR X 

9r (*>*i) v RL g^' A {ti,t 2 ) i£k (*2,*0- ( 31 ) 

Here, we have written explicitly the time convolutions 
for the sake of clarity and we have omitted the fy, in- 
tegrations, since they do not affect the time convolu- 
tions. In this expression every quantity is a 4 x 4 ma- 
trix in Nambu and spin space, and from now to the 
end of this section we get rid of the superindexes _R, A, 
since the equations for these two components are for- 
mally identical. We also drop the oo subindex, since 
all propagators in ( |3l| ) are propagators of the separated 
electrodes. The electrode Green's functions entering Eq. 
@ take the form: if) = uj(t)gj(t-t/)Uj{1f), where 

j = R,L and Uj(t) = exp [1^(^)7^3/2], <j>j(t) being the 
phase of the jth superconductor. In this expression, 
9j{t) — I ffj( e ) ex P( — de/2ir. 

We use the transformation generated by Uj (t) to trans- 
fer the time dependence from the Green's functions to the 
hopping elements 

i L R(t,t') = v LR (t)S(t - t') + J dtx J dt 2 v LR (t) 

g R (t - h) v RL (h) g L (h - i a ) i LR (t 2 ,t'), (32) 

where v LR (t) = UL(t)vLRU R (t) = v RL (t) = 
v exp [i<j>(t)T3/2]. One can easily show that all physical 
properties of the system are invariant under this gauge 



8 



transformation. Thus, we shall usually consider the T- 
matrix equation in this gauge, i.e. in which the hopping 
elements are time-dependent and the electrode Green's 
functions only depend on the time difference. 

In order to solve Eq. ( |32| ) it is more convenient to work 
in energy space where it becomes an algebraic equation. 
Thus, we Fourier transform the T-matrix with respect to 
the temporal arguments 



tLR(t,t') I jd,', 



t LR {e,e'). (33) 



It is easy to convince oneself that, due to the special 
time dependence of the coupling elements, the T-matrix 
admits a Fourier expansion of the form 

t LR (t,t') =53e^(0/a [ *Z e -fcCt-f) t LR (e,e + neV). 

n J n 

(34) 

In other words, Fourier transforming Eq. d32| ) one can 
show that tLR{e,e') satisfies the following relation 

i LR {e, = tLR(e, e + neV) S(e - e' + neV). (35) 

n 

The problem of the calculation of the current can be 
reduced to the evaluation of the Fourier components 
t n m( e ) = tLn(e+neV 1 e+meV). As it can be seen Fourier 
transforming Eq. ([32|), these components fulfill the fol- 
lowing (infinite) set of algebraic linear equations 

trim ^nm^n,m±l ~t~ ^n^nra ~t~ Vn,n— 2^n— 2,m ~t~ ^ / 'n,n+2^n-t-2,? 

(36) 

where v m -i, m = v(i + f 3 )/2, v m +i,m = v f (l - t 3 )/2, 
and the matrix coefficients £ n and V n ,m can be expressed 
in terms of the Green's functions of the uncoupled elec- 
trodes, as 

£ - ( V 9R,n+l « f 9L,n V 9R,n+l « f !h,n 
V « f 9R,n~\ V f L ,n V* 5fl : „-l v g L:7l 



V n ,n+2 = ~vfR,n+lV 



fL,n+2 9L,n+2 





9L,n-2 JL,n-2 J 



In these equations the short-hand notation g+ n = <7i(e + 
neU) for the 2x2 spin-dependent propagators has been 
used. Notice that the set of linear equations (|3^) are 
analogous to those describing a tight-binding chain with 
nearest-neighbor hopping parameters V n , n +2 and V„ jra _2. 
A solution can then be obtained by standard recursive 
techniques (see Ref. [ [iTj for details) . 

Finally, the /^.-dependence of the Green's function in 
Eq. (36) depends on our choice of the contact model. 
Thus for instance, for the disordered case the Green's 
functions appearing in Eq. ([56|) are the angle averaged 
ones, while for the case of a momentum conserving con- 
tact we must include the trajectory dependent Green's 
functions (see Eqs. 13-14). 



A. Current at finite voltage 



As commented in a previous section, the current contri- 
bution from a given trajectory may be calculated directly 
by integrating the transport equation (Q) along the tra- 
jectory over the discontinuity given by the source term. 
Thus, the time-dependent current reads as 



j(t) = eN F (j(fr,t)) t 



2V 



(38) 



where the contribution of given trajectory with momen- 
tum can be written as 



(39) 



(37) 



This expression can be greatly simplified as follows. 
First, the Keldysh components of the T-matrix can be 
eliminated in favor of the advanced and retarded com- 
ponents using the relation t K = t R ® g^ ® t A . On the 
other hand, the four elements of the enlarged space are 
not independent. For instance, it is easy to show the 
following relations: i LR = (l + t LL ® <?oo,l) ® v LR and 
i RL = v RL ® (l + goo,L ® i LL ). Using these relations it 
is rather straightforward to show that the current can be 
written ascZI 



j%, t) = Tr [f 3 (i% R ®g^,® i£ L ®gt-9L® i?R ® 9r ® 4 

+9r ® iRL ®9l® it R ~ tR L ®9l® Hr ® 9r) ] , (40) 

where we have dropped the symbol oo, since from now on in this section the only Green functions which will appear 
are the surface Green functions. 

Taking into account the Fourier expansion of the T-matrix (see Eq. [34]) , the current adopts the form 

oo 

M,t)= E SmVtie*"*®, (41) 

m= — oo 
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where the different Fourier current components can be expressed in terms of the Fourier components of the harmonics 
t nm {e) = t(e + neV, e + raeV) as follows 

Jm(i^) = / rfe^^Tr [f 3 (iLR,0n9R,n^RL,nm9L,m ~ 9L : z LR,0n9R,n^RL : nm 
J n 

+9R,0 z RL,0n9L,n l LR,nm ~ 1 RL ,0n9 L ,n z LR,nm9 R .III )] ■ (42) 

This is the general expression of the ac Josephson effect in a superconducting contact. This indicates that in 
the case of constant bias voltage there appear alternating components which oscillate not only with the Josephson 
frequency ujq, as in the case of a tunnel junction, but also with all its harmonics. When the voltage tends to zero all 
the harmonics sum up to yield the supercurrent. 

We can further simplify the expression of the current harmonics, j m . Using the equations of the T-matrix compo- 
nents it can be shown that the general relation i R '^ nm (e) = ii,']imn( € ) holds. Thus, we can simply express the current 
in terms of harmonics i^'^ nm = as follows 

jm{%) = I de^^ Tr [f 3 {t^ n g^ n i^ n gf; m — gL^on^R^mn 
J n 

,~R fAf~K fA fAf~K fA ~A \] fA n\ 
+9R,0 Z n0 9L,n T nm z n0 9 L,n z nm9 R,m ] ■ 



In order to illustrate this current formula we now inves- 
tigate the dc component for a spin singlet superconductor 
with no spin active interface in two limiting cases: tunnel 
regime and perfect transmission. For the description of a 
poorly transmissive barrier one can solve perturbatively 
the equation for t R ^ . Up to first order in the coupling 



element: t 



R,A 
LR.nm 



(V/2)(1 + T 3 )6,, 



m,n+l 



(«/2)(l 



T 3 )8 m _ n ^i and higher harmonics can be neglected. Sub- 
stituting this perturbative solution into the expression of 
the dc current, jo, one finds the traditional result 



/oo 
de p L (e- eV)p R (e) x 
-OO 

{f FD (e-eV)-f FD (e)}, 



(44) 



where f FF) is the Fermi distribution function and Pl,r 
are the local densities of states at the corresponding side 
of the interface: pi(e) — ^Im {gf(e)}. 

In the case of a perfect transmissive contact (V = 1), 
the absence of backscattering greatly simplify the calcu- 
lation of the different harmonics of the T-matrix compo- 
nents (see Ref . [ ) . In the case of a symmetric contact 
the-.de current for perfect transmission can be written 
afi 



/OO v 
de 

-°° m=0 



(I-I70IW1I 2 ) x 



{f FD (e)- f FD (e + (m+l)eV)}, 



(45) 



where 7^ = 7(5 + jeV) and 7(e) is the amplitude of an 
Andreev reflection in a N-S interface with perfect trans- 
mission. This amplitude is defined as j R = f R /{g R — i^), 
where the Green's functions are evaluated at the inter- 
face. This is a very appealing formula that tell us that 



the probability of a multiple Andreev reflection is basi- 
cally the product of the different Andreev reflections that 
take place at each side of the interface. 

In order to illustrate our approach in the case of a volt- 
age biased contact, we consider here a junction between 
two d-wave superconductors. Let us analyze in partic- 
ular the symmetric junction mentioned in section IIIB, 
whose description is based on the Green's functions of 
Eq. (20). Again, we neglect pair breaking effect and as- 
sume constant order parameter up to the interface. The 
proper self-consistent treatment of these junctions will 
be the subject of a forthcoming publication. In this con- 
tact geometry the zero-energy bound states in each side 
of the interface strongly affect the transport through this 
system. Of course, the results for the current depend on 
the contact model used. Let us first consider the case 
of a disordered contact. In this case, the propagators 
which enter in the current formula are the angle aver- 
aged ones. This implies that the anomalous propaga- 
tors vanish, which means that the current is only due to 
single-quasiparticle processes. Thus, the current formula 
reduces to 



Jo(V) = 



/OO 
de T (e, V) [f FD (e - eV) - f FD (e)] , (46) 
- 00 



where T(e, V) is an energy and voltage dependent trans- 
mission coefficient given by 



^\v\ 2 {PL(e-eV)) h { P R{e)) 



\l-\v\Hg L {e-eV)) % (g R {e))^r 



(47) 



where (pi(e))j^ (i = L,R) is the local density of states 
at the interface. 

In Fig. U we show the current-voltage characteris- 
tics and differential conductance for different values of 
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the normal transmission coefficient T> for this disordered 
model. The only abrupt feature exhibited by the cur- 
rent inside the gap occurs in the tunnel regime. The 
resonant tunneling through the zero-energy bound states 
leads to a zero-bias anomaly and the subsequent nega- 
tive differential conductance. As it is well known, the 
position and the height of the peak in the conduct ariec 
depends on intrinsic width of zero-energy statesE3L3~E3. 
It is known theoretically that the elastic scattering with 
bulk impuritiesEU or a diffusive surface layerE3 provide an 
intrinsic broadening, which for the case of Born scatteres 
is cx VTA, where V — l/2r is the effective pair-breaking 
parameter locally at the surface. In Fig. || we have in- 
troduced a small phenomenological broadening of 10~ 2 A 
to mimic this intrinsic effect. 




FIG. 3. Zero temperature I-V characteristics (upper panel) 
and differential conductance for different transmissions (lower 
panel) of a disordered contact between two d-wave supercon- 
ductors. The misorientations are a = n/4. The conductance 
is normalized by the normal state conductance and the volt- 
age is expressed in units of the maximum gap A. 



r|=0.01A 



r|=0.1A 



With this simple calculation we can already point out 
the following. Many calculations of the I-V curves in 
junctions of unconventional superconductors are limited 
to the tunneling regime, and make use of the traditional 
tunnel formula (see Eq. [S|). However, if bound states or 
divergencies are present in the density of states, this per- 
turbative expression does not even give the correct result 
for a low transparent contact, and one needs to include 
high order processes as in Eq. ([h]). Indeed, this is a 
well known problem in the context of s-wave supercon- 
ductors (for a discussion of this problem see for instance 
Ref. [ |47]]). In order to illustrate this fact, in Fig. [| we 
compare the result of Eq. (46) and the tunneling result 
for a small transparency V = 0.01. One can see in Fig. 
^(a) that there is a clear difference between these two 
results. In particular, in the low voltage regime the com- 
plete expression renormalizes the unphysical divergencies 
obtained with the tunnel formula. As the broadening of 
the level is increased the difference between these two 
results diminishes progressively (see Fig. ||(b)). Inc 
recent experiments in afe-plane tunneling into YBCOE! 
suggest that the broadening of the zero bias anomaly can 
be as large as 25% of the gap A (see Ref. [ In this 
case, there is no significant difference in the resulting cur- 
rent between both expressions. 
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eV/A 
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FIG. 4. Comparison between the tunnel formula (Eq. 
fli] ) and the exact result of Eq. (^) for a transmission 
T> — 0.01 and two different values of the broadening 77: (a) 
77 — 0.01A, and (b) 77 = 0.1A. The upper panels show the 
zero-temperature I-V curves and the lower ones the corre- 
sponding differential conductances. 
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FIG. 5. (a) Trajectory resolved zero-temperature I-V 
curves for the d-wave contact with a momentum conserving 
interface. The different curves correspond to different values 
of normal transmission coefficient. In order to see these curves 
in the same scale the current is normalized by the normal 
state conductance Gjv which includes the transmission coef- 
ficient. Moreover, the current and the voltage are normalized 
in units of the momentum-dependent gap. (b) Differential 
conductance for a transmission T> — 0.6. The vertical lines 
indicate the positions eV n = A p /n, n — 2, 10, as a guide 
for the eye. 



Let us now consider the case of a momentum conserv- 
ing contact, which is the usual assumption in the Zait- 
sev boundary conditions. To gain some insight into the 
final result, in Fig. || we show the contribution of an 
individual trajectory fy. The current and voltage are 
normalized in units of the gap seen by this trajectory. 
As can be observed, the current exhibits a pronounced 
subgap structure at voltages eV — A p /n, where n is a 
integer number, together with the appearance of negative 
differential conductance (this can be seen better in the 
lower panel of this figure). These features are a simple 
consequence of the resonant tunneling across the zero- 
energy bound states. Indeed, this type of TV has been 
previously obtained in the context of a junction between 
two conventional superconductors coupled by means of 
a resonant transmission (see Ref. [ |76|]77| ]). Notice also 
the presence of a zero bias peak, specially clear for low 
transparencies, and which is a consequence of the small 
broadening introduced in the calculation. The total cur- 
rent is obtained by averaging over the different trajec- 
tories. Thus, the final result depends on the model for 
the angular dependence of the normal transmission co- 
efficient. With any reasonable reasonable model most of 



the features of the trajectory resolved current disappear. 
In particular, the subharmonic gap structure is washed 
out, and it only remains a peak in the conductance at 
eV f» A (see for instance Ref. [ [HI). 



V. CURRENT FLUCTUATIONS 

During the last years it has become progressively 
clearer that a deep understanding of the electronic trans- 
port in mesoscopic systems requires the analysis of quan- 
tities which goes beyond the straightforward measure- 
ment of the current-voltage characteristics. In this sense 
the noise or time-dependent current fluctuations has 
emerged as a very useful tool which provides new infor- 
mation on the time correlations of the current, informa- 
tion about channel distributions, statistics and charge of 
the carriersES. In the case of superconducting contacts, 
most of the activity has-.been restricted to the case of 
s-wave superconductorsE2TE3. In the case of unconven- 
tional superconductors there are only a few theoretical 
works in the context of hybrid, [Structures like normal 
metal/d-wave superconductorsSE.il. We believe that in 
the next future the measurement of current fluctuations 
will be an important tool for a deeper understanding of 
the symmetry of the order parameter and origin of the su- 
perconductivity in general in the case of HTS materials. 
For this reason in this section we describe the calculation 
of the noise spectrum within our approach, valid for any 
type of contact between unconventional superconductors. 

Let us remember that the noise is characterized by its 
spectral density or power spectrum S(u>), which is sim- 
ply the Fourier transform at frequency u) of the current- 
current correlation function 



sh = / d(f - 1) e^( f '~ f )(5}(t>)5m + smm')) 

d(t'-t) eM'-O K(t,t'), 



(48) 

where Sj(t) — < ](t) > are the fluctuations in the 

current. 

In order to obtain the expression of the current-current 
correlator, first we need an expression for the current op- 
erator. Within our model this operator evaluated at the 
interface can written as follows 

= ley] |«LH,gC^ g (f)Cfl |g (t) - 

(49) 

This expression is a simple consequence of the continuity 
equation for the currenttj. 

In order to calculate the noise we need in principle 
to evaluate correlators of four field operators. However, 
we are working in the framework of a mean field the- 
ory, which means that we can decouple these correlators 
in terms of one-particle Green's functions using the Wick 
theorem. With this in mind, it is straightforward to show 
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that the kernel K(t,t') can be expressed in tepp 
interface Keldysh Green's functions as followsEil 

K(t,t') = e 2 {TV [v RL G< L (t,t')v LR G> R (t',t) 
+v LR G< R (t,t')v RL G> L (t',t) 

-VR.LGf (t,t')v RL G> R (t',t) 



is of the 



trace Tr and the matrix T3 which act in the "reservoir" 
space. Then, the noise kernel reads 



K(t,t') 



-e 2 fr 



vG < {t,t')nvG > {t',t)+ 
vG > (t,t , )f 3 iG < (t , ,t) 



(52) 



-«LflG< x (t, 0«l«<% l (*', t)] + (t - 0} . (50) 

where the functions G < and G > are related to the usual 
advanced, retarded and Keldysh functions in the follow- 
in cr wmr 



G< = (G K - G R + G A ) /2 



(51) 



G> = (G x + G fl - G A ) /2. 
In order to compactify the notation, we introduce the 

This expression allows us to write the noise kernel as follows 



Now, in order to eliminate the Green's functions in 
favor of the T-matrix elements, we use the relation 

G<'> = f 1 + G R o vj o G<-> o (l + v o G A ) , (53) 

GA{e)-G*{e)]fMe) and G>(e) = 
(/fo(e) — 1). Now, making use of Eqs. 
(3-7) it is easy to see that the following relation holds 

ioG<<>=f R o G ^>o(l+f A oGi)). (54) 



where G^(e) = 
G A (e)-G R (e) 



A-(t,0 



VTr • 



T K oG< o 1+T^oG^ 



(*, 7* 



T o G^, o(l+ T o G 



)](«',*) + (*-*')}• (55) 



As explained in section II, once we have eliminated the full Green's functions in the noise kernel, we can perform 
the quasiclassical ^-integration and substitute the quasiclassical Green's functions, goo, for the full ones, Goo. Thus, 
the noise kernel can be expressed finally as 

K(t,t')=-e 2 fr{[f R ®g< ® (l + f A ® gi)] (t,t') f 3 [f R <g> g> ® (l + f A ® ~g£)] (t' ,t) + (t -+ *')} . (56) 

Let us stick now to the case of a constant bias voltage applied across the interface. In this case, for both contact 
models considered in section III, we can resolve the current fluctuations in trajectories as follows 



S(w,t) = e 2 N F (S(&,u;,t))^, 
where the time-dependent contribution of given trajectory with momentum can be written as 



(57) 



(58) 



m— — 00 



where the different ac-components of the noise can be expressed in terms of the Fourier component of the T-matrix 
elements in the following way 



dc 



EM 



f R ~g< [1 + T*g, 



A~A 

DC 



On 



T 



R~> 



(l + f A g£)] (e + uj) + (e^e + Lo)\. (59) 

\ /J nm J 



Notice that in the case of a junction between two super- 
conductors, the noise, as the current, oscillates on time 
with all the harmonics of the Josephson frequency. No- 
tice also that we have reduced the calculation of this 
quantity to the determination of the different Fourier 
components of the T-matrix elements, which has been 
detailed in section IV. 

In order to illustrate the calculation of the current fluc- 



tuations, we consider the contact between d-wave super- 
conductors analyzed in the previous section. In particu- 
lar, we present results for the zero-frequency noise at zero 
temperature, S, i.e. the zero-frequency shot noise. At 
this point, it is worth remarking that by zero- frequency 
noise we mean noise at a frequency lower than any rele- 
vant energy scale in our problem—gap for instance, and 
high enough to neglect 1// noiseo. Again, the final re- 
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suit depends on the type of contact model under inves- 
tigation. Let us start discussing the shot noise for the 
disordered contact. In this case, due to the vanishing of 
the anomalous Green's functions, the whole calculation 
reduces to the determination of the quasiparticle contri- 
bution, which in term of the transmission coefficient of 
Eq. ( f47| ) can be written as 



S 



eV 



deT(e,V)(l-T(e,V)). 



(60) 



This is simply the result that one obtains for a normal 
contact with aa-energy and voltage dependent transmis- 
sion coefficient^! In Fig. ^| we show the result of Eq. 
( |60| ) for different normal transmissions. In the tunneling 
regime the shot noise is S(V) « 2ej(V) and the most 
remarkable feature is the zero bias anomaly (see inset 
of Fig. H). In the case of perfect transmission there 
is a non-zero noise due to the fact that the transmis- 
sion coefficient T(e, V) is less than one in the gap region. 
For voltages much greater than the maximum gap, then 
T(e ^S> A,V = 1.0) — * 1, what makes that the noise at 
V = 1.0 saturates in the high voltage regime. 




FIG. 6. Zero-frequency shot noise for the disordered con- 
tact considered in Fig. H. The inset show the low bias limit 
of the curve in the tunneling regime (T> — 0.01). 

More interesting is the case of the momentum conserv- 
ing interface. In Fig. |we show the contribution of a tra- 
jectory of momentum f^,. As in the case of the current, 
the shot noise exhibits a rich subharmonic gap structure, 
which persists almost up to perfect transmission. The 
shape of the different curves can be understood in the 
same terms as the BCS case (see Ref. [ |8l|), with the ad- 
ditional ingredient of the resonant tunneling through the 
zero energy states. Of course, as in the case of the cur- 
rent, most of these features disappear after performing 
the angular average. 
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FIG. 7. The upper panel shows the trajectory resolved 
zero-frequency shot noise for the momentum conserving case 
considered in Fig. ^. The shot noise and voltage are normal- 
ized by the trajectory-dependent gap A p . In the lower panel 
one can see the effective charge defined as Q* = S/2ej as a 
function of voltage for a transmission T> — 0.6. The verti- 
cal lines indicate the position of the voltages eV„ = A p /n, 
n = 2, 10. 

In the case of conventional superconductors, the shot 
noise has been proposed as a tool for measuring the 
multiple charge-quanta transferred by the multiple An- 
drev reflections^. Obviously we can pose here the same 
question in the case of unconventional superconductors. 
Indeed, a noise experiment—has been recently proposed 
by Auerbach and Altmar£3 to discriminate between 
two possible explanations of the pronounced subhat- 
monic gap structure observed in YBCO edge junctionsEij. 
Namely, usual multiple Andreev reflections in a d-wave 
superconductor and magnon pair creation in the context 
of the SO(5) theory. In this latter case the observed 
charge should be Q* = 2ne, where n — 1,2, at a volt- 
age eV^ = A/n. This result has to be compared with 
Q* = ne expected in the traditional view of MAR. In or- 
der to contribute to the solution of this puzzle, we show 
in Fig. (lower panel) the effective charge, Q* — S/2ej, 
for a transmission T> = 0.6. This result confirm the tradi- 
tional interpretation that in the MAR process of order n 
a charge ne is transferred. Usually, in order to observe a 
clear quantization of the charge one should go to the tun- 
neling regimeEll, but in this case this is not necessary due 
to the resonant tunneling through the zero-energy states. 
Notice again that this is the contribution of a single tra- 
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jectory and after angle averaging this clear quantization 
of the charge with voltage disappears. The exhaustive 
analysis of the shot noise in d-wave contacts will be pre- 
sented in a forthcoming publication. 



VI. CONCLUSIONS 

We have shown how a Hamiltonian approach and the 
quasiclassical theory of superconductivity can be com- 
bined to give a powerful tool to analyze electronic and 
transport properties of superconducting junctions. In 
particular, we have demonstrated that a simple Hamilto- 
nian description of an interface can be used to model a 
great variety of contacts. This Hamiltonian description 
can be brought into quasiclassical theory via a T-matrix 
equation, resulting in a new formulation of boundary con- 
ditions. These boundary conditions do not contain any 
spurious solution and can be efficiently solved to com- 
pute any transport property. The broad applicability of 
this formulation covers cases ranging from conventional 
superconductors to unconventional ones, clean systems 
and diffusive ones. Moreover, it can be applied to spin 
active interfaces and it is well suited for the description 
of time-dependent phenomena like the I-V characteris- 
tics and the noise properties of junctions with arbitrary 
transmission and bias voltage. We have illustrated this 
approach with the calculation of Josephson current in a 
great variety of situations. The calculation of I-V char- 
acteristics and the noise has been exemplified with the 
analysis of a contact between two d-wave superconduc- 
tors. In particular, we have briefly discussed the use of 
shot noise as a possible tool for measuring the charge of 
the Andreev reflections in unconventional superconduc- 
tors. 
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